Anomalous spin susceptibility and magnetic polaron formation in the double exchange 

systems 
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The magnetic susceptibility and spin-spin correlation of the double-exchange model for doped man- 
ganites are investigated through the Monte Carlo calculations on the three-dimensional lattice model. 
Deviations of the susceptibility from the Curie- Weiss behavior above the ferromagnetic ordering tem- 
perature T c seem to indicate a formation of local ferromagnetic clusters in the vicinity of T c , which is 
consistent with recent electron paramagnetic resonance experiments for La 2 /3Cai/3Mn03. A further 
analysis of the spin-spin correlations show the ferromagnetic cluster size to be three-to-four lattice 
spacings, suggesting that the charge carriers may form magnetic polarons. 



PACS numbers: 75.30.Et, 75.30.-m, 71.10.-w 
I. INTRODUCTION 

Recent interests in mixed- valence manganites with the 
chemical compositions of Ri- x A^MnO^ (where i?=rare 
earth; A=Ga,, Sr, Ba) are largely ascribed to their poten- 
tial technological applications E The most notable fea- 
ture of these materials is an extremely large change of 
resistivity under the application of a magnetic field near 
ferromagnetic ordering temperature T c , which is known 
as colossal magnetoresistance (CMR). A metal-insulator 
(MI) transition accompanied by the ferromagnetic- 
paramagnetic phase transition is occurred near T c . The 
strong connection between the MI transition and the fer- 
romagnetic spin alignment has been understood in terms 
of the double-exchange (DE) mechanism. Era The conduc- 
tion electrons in the e g orbital of Mn 3+ ions are hop- 
ping in the background of Mn 4+ {t^g) ion spins with 
an experience of a strong on-site Hund's rule coupling, 
thereby leading to an effective hopping integral of the 
form tcos(9ij/2) where Oij is the relative angle between 
Mn 4+ ions. 

Although the electron-lattice interaction arising from 
the dynamic Jahn- Teller distortion is considered to be 
important for the understanding of overall trends of CMR 
phenomenau, the lattice polaron formation is incomplete 
to explain the transport properties in connection with the 
observed CMR phenomenaJjLl Contrary to the above ap- 
proach, it is suggested that the e g carriers can be trapped 
by spin-disorder scattering due to the local deviations of 
the ferromagnetic surroundings, resulting in thx: forma- 
tion of magnetic polarons ia the vicinity of T c .ErEl Indeed, 
a lot of theoretical workstjli 4 ! have demonstrated that 
the magnetic polaron formation plays a crucial role in 
CMR phenomena in the paramagnetic state. 

From the experimental point of view, much efforts have 
been devoted to understand their unique magnetic fea- 



tures in the paramagnetic state. For instance, the high 
temperature inverse susceptibility for Lao.7Cao.3Mn03 is 
smaller than expected frop-J-he Curie- Weiss law at low 
magnetic field H < 0.1 T.EjI1j The deviation of suscep- 
tibility from Curie- Weiss law is also found in the layered 
material of Lai.35Sri.65Mii207, suggesting the formation 
of ferromagnetic cluster above T c .li3_Recent electron para- 
magnetic resonance (EPR) studyLa revealed that EPR 
signals decrease exponentially at temperatures slightly 
above T c , which results in an enormous enhancement of 
the effective spins of about ~ 30 per formula unit —More- 
over, small angle neutron scattering experimentE3 esti- 
mated the size of the ferromagnetic cluster in the param- 
agnetic state to be Acp-f 12A. There has also been an in- 
creasing realizationE3 c3 which manifest the importance 
of the magnetic fluctuations that is beyond the mean- 
field prediction in the paramagnetic state of the doped 
perovskite manganites. 

The main purpose of this paper is to elucidate the 
characteristic features of the ferromagnetic polaron in the 
vicinity of T c . In order to account for the spin fluctuati 
effects accurately, we adopt a Monte Carlo methocO 
for the DE model on the three-dimensional (3D) lattice. 
We found that the temperature dependence of inverse 
susceptibility deviates from the expected Curie- Weiss be- 
havior, which is in good agreement with experimental 
observations. Also, the activation energy estimated from 
the non- Curie- Weiss part of the susceptibility coincides 
with that of EPR measurements for La2/3Cai/3Mn03.t3 
From these results together with the calculated spin-spin 
correlation, we suggest that charge carriers above T c form 
the ferronragnetic cluster with the size of three-to-four 
lattice spacings. 



II. MODEL AND CALCULATIONS 

One of the simplest models for the description of the 
paramagnetic-ferromagnetic phase transition in doped 
manganites is a single-orbital DE model Hamiltonian. In 
the strong-,Iiund's coupling limit (J# — > oo), it can be 
written asETt3 



" = -£( 



tijcfci + h.i 



-/*£#* 



(i) 



where the operator c[ creates a spinless conduction elec- 
tron at site Ri and h is an external magnetic field. The 
hopping amplitude in the strong Hund's coupling limit is 
given by 



tn = t[ COS — COS-r- 



(2) 



where the localized tig spin Si is treated as a classi- 
cal spin Si = |iSi|(sm#j cos fax + sin #j sin fay + cosOiz) 
with the polar angles {9i,fa} characterizing the orien- 
tation of the localized spin Si since the quantum ef- 
fects can be neglected near room temperature. Further- 
more, we assume that the azimuthal angle fa rotates in- 
dependently since the contribution of the phase factor 
in the hopriLag integral to the partition function can be 
neglected. EjEj Thus, it is possible to replace ty by the 
DE form tij = icos(0y/2), where 9ij is the relative angle 
between neighboring Mn ions, and the conduction elec- 
trons are treated as a quantum object. 

The grand canonical partition function of the present 
model with chemical potential can be denoted by 



Z = TrcTr f 



,-0(H-iiM) 



(3) 



where Tr$ and Tr c represent traces over the localized 
spin configuration and the conduction electron degree 
of freedom, respectively, and TV = J2i c l c i- Here the 
conduction electron density (n) — {J\f)/L 3 is determined 
by adjusting chemical potential /i. Since the DE model 
Hamiltonian is quadratic with respect to the Fermion 
operator {ci}, Tr c is directly calculated by diagonalizing 
the electronic part of the Hamiltonian. By performing 
the diagonalization of the (L 3 x L 3 ) hermitian matrix 
for each given spin configuration {9i,fa}, we obtain the 
L 3 eigenvalues denoted by e a (8i,fa). Thus the resulting 
partition function becomes 



Z = X\( dfliSinfli / dfa) TT(1 
i Vio Jo J a=1 



+ e -/J(ea-*.)). 



(4) 



Now, we can apply a Monte Carlo integration procedure 
for the summation over the configuration angles {9i, fa} 
of localized spins using a standard Metropolis algorithm. 



The actual calculations are performed for 3D cubic lat- 
tices L 3 = 6 3 with periodic boundary conditions in spa- 
tial directions. Typically, we take 5000 Monte Carlo 
steps per site for statistical average. Thermodynamic 
quantities of interest are obtained directly from the ther- 
mal average of spin configurations and the eigenvalues of 
Hamiltonian. The carrier density (n) « 0.5, i.e., the hole 
density x = 1 — (n), obtained by fixing ^ = 0.0. 



III. RESULTS 

Before studying the magnetic fluctuation effects, we 
first investigate the nature of the magnetic transition in 
the ferromagnetic state. Figure 1 shows the temperature 
(T) dependence of the magnetization, M = (J^\ Sf )/L 3 , 
denoted by solid squares. As clearly shown in Fig. 1, 
the magnetization drops sharply when the ferromagnetic 
ordering temperature reaches to T c « 0.125i and the 
so-called ferromagnetic-to-paramagnetic transition takes 
places near T c . The dashed line represents the mean-field 
prediction of DE model by Kubo and Ohata, leading to 
M = M \T C g- T] 1 / 2 with M = ^/b/Tc for S = 2 as T c is 
approached. U This magnetic transition is closely related 
to the other transition of apparently different character, 
i.e., the T-dependence of resistivityu The obtained Monte 
Carlo result, however, varies more slowly than that of the 
mean-field prediction of DE model. Moreover it is noted 
that the magnetization data can be fitted with a power 
law of M(T) - \T C - Tf just below T c . The solid line 
in the Fig. 1 is a fitting curve, yielding j3 — 0.32 ± 0.01. 
This value of the exponent iLis very close to that of 
Lao.7Sro.3Mn03 single crystalE 3 ] and 3D Heisenberg class 
which is j3 = 0.33 in the vicinity of T c . Above T c , how- 
ever, our calculated result of M is different from zero 
due to the finite-size effects. The similar results are also 
obtained by using different Monte Carlo methods for a 
larger size of unit cells 20 x 20 x 20t3 
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FIG. 1. Temperature dependence of magnetization M(T). 
The dashed and the solid line is the mean-field prediction of 
double-exchange model and a fitting curve proportional to 
\T C -Tj ' 32 , respectively. 

Now we turn to study the magnetic fluctuation effects 
in the paramagnetic state. The magnetic susceptibility 
directly measures the magnetic fluctuation as shown in 
the following equation: 



dM 
~~dh 



() 



Tr y\ gz^HCH-^M-hY,^!) 



L 3 dh TTe -i3(H-p,Af-hY; z S*) 
(M 2 ) - (M) 2 



T 



(5) 



Figure 2 provides % as a function of T/T c denoted by 
open squares. With increasing temperature, \ increases 
to a maximum near T c and then quickly decreases above 
T c . It should be noted that the susceptibility is well de- 
scribed by a power law of \ ~ \T— T c |~ 7 as T approaches 
close to the critical temperature T c in the paramagnetic 
state, where the critical fluctuations are dominant. The 
inset of Fig. 2 shows ln(x) vs. ln(|T — T c \) curve, yield- 
ing the exponent 7 = 0.92 ± 0.01 and the dashed line 
in the Fig. 2 is a fitting curve. This value of the critical 
exponent 7 is smaller than those ofihe experimental sug- 
gestion of 1.22 for Lao.7Sro.3Mn03c3 and-the mean-field 
theory of Curie- Weiss law which is 7 = lc3. However, it 
is not unreasonable if we take account of the fact that 
the present value is determined only from the 6 3 finite 
size cluster calculation with no scaling analysis. 
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FIG. 2. Temperature dependence of magnetic susceptibil- 
ity x- The dashed line is the fitting curve of \T — T c | -7 with 
7 = 0.92. The vertical lines indicate the estimation of error 
bar. Inset shows ln(x) vs. ln(|T — T c \) above T c . 



d(x~ 1 )/dT in the high temperature regimecS of a typical 
itinerant ferromagnetic such as nickel (T c = 623 K) is 
known to be temperature-dependent. In order to under- 
stand the spin entities in the paramagnetic state, we plot 
1/X vs. T/T c in Fig. 3. In the high-temperature region of 
T > 2.5T C , the corresponding l/x follows a linear Curie- 
Weiss behavior where each spin has a non-interacting 
magnetic moment. The solid line is the fitting curve of 
the Curie- Weiss law, \ ~ C/(T — 8), yielding the mean- 
field Curie temperature 6 = 0.5T C and C = 0.0008. How- 
ever, starting from and below 2T C , l/x shows a distinct 
deviation from the expected Curie- Weiss law, suggesting 
a possible presence of short-range ferromagnetic clusters 
due to the spin fluctuations above T c . This feature are 
consistent with observations in the recent experimental 
measurements on the single crystal of Lap.yCao.sMnOsEJ, 
the thin film of Lao.6Yo.o7Cao.3*Mn03E3, and the layered 
material of Lai.35Sri.65Mn207tZl in low magnetic field 
h < 0.1 T. For the high magnetic field, however, l/x be- 
comes larger than the one expected from the Curie- Weiss 
lawE£rE£|, indicating a suppression of ferromagnetic fluc- 
tuation due to the carrier delocaliaartion and magnetic 
ordering by the applied magneticJ13'E2l 

Assuming that ferromagnetic clusters exist, one should 
see an effective moment which is larger than that of the 
appropriate average of Mn 3+ and Mn 4+ moments. From 
the consideration that P neighboring sites with ferromag- 
netically aligned spins form a spin polaron of the spin 
(Si + PS2) due to the double exchange mechanism, the 
effective spin for the fluctuation effects becomesQ 

S 2 ff = x(S 1 +PS 2 ){S 1 +PS 2 + l) + (l-x-Px)S 2 (S 2 + l), 

where Si and S 2 identify the spin of Mn 4+ and Mn 3+ 
species, respectively. This estimation of the square of 
total spin will become larger than the one without spin 
polaron formation, i.e., P — case, and consequently 
the inverse susceptibility will be smaller compared to the 
Curie- Weiss law. These features are consistent with our 
current Monte Carlo results as,_well as other results of 
the exact diagonalization studyEI in which the effective 
spins of e g carriers gpow from 1/2 to 7 and the EPR 
experiments above T C .EZI Thus, it is reasonable to inter- 
pret that the spin polarized carriers form ferromagnetic 
polarons with each individual cluster retaining a large 
magnetic moment. 



On the other hand, in the high temperature regime of 
~ 2T C , the magnetic susceptibility of the DE system ex- 
hibits a peculiar temperature dependence. In fact, these 
magnetic entities in the paramagnetic state have been 
investigated over a long period of time. For instance, 
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FIG. 3. Temperature dependence of the inverse suscepti- 
bility 1/x- The solid line is the best fit of the Curie- Weiss 
behavior of the form x = C/(T - 0.5T C ). 

It is interesting to compare our Monte Carlo results 
with the outcomes from the EPR measurements. For 
this purpose we define the deviation of susceptibility from 
the Curie- Weiss law as Ax = x — C/(T — 9). In the 
intermediate temperature region of T c < T < 2.5T C , we 
find that Ax can be fit by the following form 



A X = Xoexp(£ a /T) 



(G) 



where E a is an activation energy which scales to the T c . 
Figure 4 shows ln(Ax) vs. 1000/T for fixing T c = 270 K 
denoted open squares. The solid line is a fitting curve, 
yielding E a = 10T C which corresponds to 0.22 eV for 
t = 0.17 eV, when Xo = 9 x 10~ 6 is used. 

It is interesting to note that the EPR inten- 
sity decreases exponentially above T c and has a 
strong correlation with the deviation of the magnetic 
susceptibilityOtJ. For a comparison with our result of 
Ax, the log of intensity as a function of 1000/T for 
La 2 /3Ca 1/3 Mn03 (T c w 270 K) from Ref. [18] is shown 
in the inset of the Fig. 4. The shape of the T-dependence 
of EPR intensity is quite similar to the Monte Carlo re- 
sult of Ax- The E a obtained from EPR measurement is 
approximately 7T C which is slightly smaller than that ob- 
tained-from the Ax- According to the studies of Oseroff 
et al.f3 those spin entities are associated with some form 
of the magnetic cluster with the spins of about ~ 30 per 
formula unit. In other words, the E a would be required in 
order to dissociate these spin entities made of collection 
of individual spins. Therefore, it is suggested that the 
magnetic polaron formation driven by the spin-disorder 
scattering should be generic of the MI transition regimes 
and is likely associated with the CMR phenomena in the 
doped manganites. 
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FIG. 4. ln(Ax) vs. 1000/T for T « 270 K and the solid 
line represents the best fit to Eq. (6). In the inset, EPR inten- 



sity of ln(7) vs. 
from Ref. [18]. 



1000/T for La2/3Cai /3 Mn03 (T c as 270 K) 



Finally, we estimate the effective sizes of the ferromag- 
netic clusters in the paramagnetic state. Figure 5 shows 
the spin-spin correlation (JT Si ■ Si +r ) as a function of 
T/T c in the vicinity of T c for several values of lattice 
sites. The shortest correlation for r = 1 is robust even 
for T > 2T C and, in particular, almost 15% of the maxi- 
mum at T c . Furthermore, the short-range ferromagnetic 
correlation is clearly seen in the intermediated regime of 
T c < T < 2T C . The size of the ferromagnetic surround- 
ings is of the order of 3 ~ 4 lattice spacing distances, i.e., 
12A^16A if we consider a size of lattice as 4A, which, 
is in good agreement with recent neutron scatteringHj'Ell 
Moreover, a simple estimation of the ferromagnetic corre- 
lation length, P fa (oi/T) 3 / 5 , in terms of a ferromagnetic 
spin-polaron picture has been introduced by VarmaB. At 
the temperature range of paramagnetic state, in partic- 
ular at room temperature ~ 300 K, the size of the spin 
polaron is estimated to be a few lattice size for t ~ 0.2 eV, 
which agrees with the current results. It implies that the 
spin-polarized carriers can be trapped into a local fer- 
romagnetic surroundings due to the spatially fluctuating 
spin correlations, resulting in the formatioa-of the mag- 
netic polarons in the paramagnetic phase.Bli2l 
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5. Temperature dependence of the spin correlation 
• Si +r ) for x = 0.5 and L 3 = 6 3 . 



IV. CONCLUSIONS 



Using Monte Carlo methods, we have studied the mag- 
netic fluctuation effects near and above T c within the 
framework of DE model. We found that the temper- 
ature dependence of the susceptibility shows a sharp 
peak near T c , and the inverse susceptibility displays non- 
Curie- Weiss behavior above T c . The activation energy 
obtained from the deviation of susceptibility from the 
Curie- Weiss behavior is consistent with the EPR mea- 
surements. These results clearly demonstrate the forma- 
tion of the magnetic polaron with short-range ferromag- 
netic ordering in the paramagnetic state. Moreover, the 
ferromagnetic correlation length is estimated to be 3 ~ 4 
lattice spacings which is in goo d . ag reement with recent 
neutron scattering experimentsE3n9. From the results, 
it is suggested that the magnetic polaron formation is 
responsible for the magnetic transition and the magneto- 
transport properties in doped CMR manganites. 
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